Time Delay Interferometry for LISA with one arm 
dysfunctional 



S. V. Dhurandhar 1 , K. Rajesh Nayak 2 and J-Y. Vinet 3 

1 IUCAA, Postbag 4, Ganeshkind, Pune - 411 007, India. 

2 IISER-Kolkata, PO: BCKV Campus Main Office, Mohanpur - 741252, India. 

3 ARTEMIS, Observatoire de la Cote d'Azur, BP 4229, 06304 Nice, France. 



Abstract. 

In order to attain the requisite sensitivity for LISA - a joint space mission of the ESA and 
NASA- the laser frequency noise must be suppressed below the secondary noises such as the 
optical path noise, acceleration noise etc. By combining six appropriately time-delayed data 
streams containing fractional Doppler shifts - a technique called time delay interferometry 
(TDI) - the laser frequency noise may be adequately suppressed. We consider the general 
model of LISA where the armlengths vary with time, so that second generation TDI are 
relevant. However, we must envisage the possibility, that not all the optical links of LISA will 
be operating at all times, and therefore, we here consider the case of LISA operating with 
two arms only. As shown earlier in the literature, obtaining even approximate solutions of 
TDI to the general problem is very difficult. Since here only four optical links are relevant, 
the algebraic problem simplifies considerably. We are then able to exhibit a large number 
of solutions (from mathematical point of view an infinite number) and further present an 
algorithm to generate these solutions. 



PACS numbers: 95.55.Ym, 04.80.Nn, 07.60.Ly 
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1. Introduction 

LISA - Laser Interferometric Space Antenna - is a proposed ESA-NASA mission which will 
use coherent laser beams exchanged between three identical spacecraft forming a giant (almost) 
equilateral triangle of side 5 x 10 6 kilometres to observe and detect low frequency cosmic GW 

DP- 

Laser frequency noise dominates the other secondary noises, such as optical path noise, 
acceleration noise by 7 or 8 orders of magnitude, and must be removed if LISA is to achieve the 
required sensitivity of h ~ 10 -22 , where h is the metric perturbation caused by a gravitational 
wave. In LISA, six data streams arise from the exchange of laser beams between the three 
spacecraft approximately 5 million km apart. These six streams produce redundancy in the 
data which can be used to suppress the laser frequency noise by the technique called time-delay 
interferometry (TDI) in which the six data streams are combined with appropriate time-delays 
[2]. A mathematical foundation for the TDI problem for the static LISA was given in [3], where 
it was shown that the data combinations cancelling laser frequency noise formed the module of 
syzygies over the polynomial ring of time-delay operators. For the static LISA, the polynomial 
ring was in three variables and also commutative. The generators of the module were obtained 
via Grobner basis methods. This scheme can be extended in a straight forward way to include the 
(Sagnac) effect arising from the rotation of LISA, where the up-down optical links are unequal 
and as a consequence, now one has six optical links, but the armlengths are still constant in time. 
These are the modified 1st generation TDI. Here the polynomial ring is still commutative although 
over six indeterminates. The generators of the module were obtained in [4]. In the general case 
of time varying armlengths, the polynomial ring is in six variables (the six optical links) but now 
it is noncommutative. The algebraic problem in the general case has been discussed in [5], but 
the solution to the general problem seems extremely difficult. We still have a linear system which 
leads to a module (a left module), but it seems difficult to obtain its generators in general. In 
the non-commutative case, even the Grobner basis algorithm may not terminate. 

We must envisage the possibility that not all optical links of LISA will be operating at all 
times for various reasons like technical failure for instance or even the operating costs. Therefore, 
it is important to discuss the question when not all the links operate. Here we look at a specific 
situation where the data is available from only two arms or four optical links. This should not 
much affect the information that can be extracted from the data, because this is essentially a 
Michelson configuration which is known to be quite useful. The practical advantage for this 
case is that the algebraic problem simplifies considerably and therefore becomes tractable. We 
can reduce the problem to that of only one linear constraint on two polynomials, although the 
equation is still noncommutative. We show that we can generate an infinity of solutions by 
using a combinatorial algebraic approach which lists all such solutions in a systematic way. The 
solutions - that is the laser frequency noise is suppressed for these time-delayed combinations - 
are approximate in the sense that the L and L 2 terms are ignored in the calculation, where L(t) 
is the generic length of the LISA's arm or the optical link and the 'dot' denotes derivative with 
respect to time. The solutions are based on vanishing commutators (in this approximation) . We 
enumerate such commutators and for each such commutator there is a corresponding solution. 
We present an algorithm to construct these solutions. 

A geometric combinatorial approach was adopted in [5] where several solutions were 
exhibited. Our approach is algebraic where several operations are algebraic operations on strings 
of operators. The algebraic approach has the advantage of easy manipulation of data streams, 
although some geometrical insight could be at a premium. The analysis presented here can be 
used also for other future space detectors like ASTROD where again a similar situation arises 



Time Delay Interferometry for LISA with one arm dysfunctional 



3 



2. Preliminaries and notation 



In this section we recapitulate the results obtained in the earlier literature already mentioned 
and also set up notation. In the literature on TDI of LISA there are several notations - here we 
choose the one which seems simplest for our purpose. We follow the notation and conventions 
of [3] and [4]. The six links are denoted by U l ,V l ,i = 1,2,3. The time-delay operator for the 
link U 2 from S/CltoS/C2orl — > 2 is denoted by x in [4] and so on in a cyclic fashion. The 
delay operators in the other sense are denoted by l,m,n; the link —V 1 from 2 — > 1 by / and 
similarly the links V 2 ,V 3 are defined through cyclic permutation. Figure [1] depicts the optical 
links as described. 



SCI 




SC3 



SC2 



Figure 1. The beams and the corresponding time-delays are shown schematically in the 
figure. The functional arms are depicted with a continuous line while the dysfunctional arm 
is shown with a dashed line. 



Let Ci(t) — Ai/i(t)/i/Q represent the laser frequency noise in S/C i. Let j be the delay 
operator corresponding to the variable armlcngth Lj(t), i.e. jCi(t) = d(t — Lj(t)). Then we 
have, 

U 1 = Ci - zC 3 , 

V 1 = IC 2 - C x . (1) 

The other links in terms of Ci(t) are obtained by cyclic permutations. Also in the U % we 
have not included contributions from the secondary noises, gravitational wave signal etc. since 
here our aim is to deal with laser frequency noise only. Any observable X is written as: 

X= Pl V l + qi U\ (2) 

where Pi,qi,i — 1,2,3 are polynomials in the variables x,y, z,l,m,n. Thus X is specified by 
giving the six tuple polynomial vector (pi,qi). Writing out the (V l ,U l ) in terms of the laser 
noises Cj(t), and in order that the laser frequency noise cancel for arbitrary functions Cj(t), the 
polynomials (pi,qi) must satisfy the equations: 

Pi - <?i + Q2X -ps,n = 0, 
Pi - q% + q3V - Pil =0, 

P3 - 93 + qiz -P2m = , (3) 
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where one must remember that the order of the operators is important when the armlengths 
are time-varying. The solutions (if they exist) to Eq. Q are the second generation TDI. One 
finds that at least in the approximation we consider they do exist. For a reasonably optimised 
model of LISA, say given in [8], L ~ 10~ 6 metres/sec 2 and thus even if one considers say 20 
successive optical paths, that is, about At ~ 330 seconds of light travel time, At 2 L ~ 0.1 metre. 
This is well below few metres and thus can be neglected in the residual laser noise computation. 
Moreover, L 2 terms (and higher order) can be dropped since they are of the order of < 10~ 15 
(they come with a factor 1/c 2 ) which is much smaller than 1 part in 10 s , which is the level at 
which the laser frequency noise must be cancelled. Thus, we keep terms only to the first degree 
in L and also neglect higher time derivative terms in L. 

The solutions to these equations are the polynomial vectors X = (j>i,qi) and they form a 
left module M. over the non-commutative ring K, = Q{x,y 1 z,l,m,n), where Q is the field of 
rational numbers and x,y, z,l,m,n are the indeterminates. A left module means that one can 
multiply a solution (pj, qi) from the left by any polynomial in 1C, then it is also a solution to the 
equations Q and therefore in the module M.. It is difficult to solve these equations in general 
for the non-commutative case. Previously we have shown in [3] that the ring K, can be reduced 
to the quotient ring K, which is defined by quotienting K, by the ideal IA generated by vanishing 
commutators. Consider commutators of the form \x\x 2 ...x n ^y\y 2 ...y^\ where n > 2 and Xk or 
y m represents any one of the delay operators x, y, z, I, m, n. Upto the order of approximation we 
work in: 

n n n n 

[x!x 2 ...x n , yiy-2-Vn] = Y L ** ^y,„ - Y L v™ Y ' ( 4 ) 

fe=l m=l m=l fc=l 

From this equation it immediately follows that if the set of operators yi, y 2 , y n is a permutation 
of the operators x\, X2, x n , then the commutator: 

[x 1 X2...x n ,yiy2...y n ]=Q, (5) 

upto the order of the approximation. The ideal IA is generated by such commutators. 

Although we may be able to apply non-commutative Grobner basis methods, the general 
solution seems quite difficult. However, we could consider a simpler problem of one arm of LISA 
being nonfunctional. Then we have a simpler situation which turns out to be tractable. 



3. Solutions with one arm nonfunctional 



We arbitrarily choose the non-functional arm to be the one connecting S/C 2 and S/C 3. This 
means from our labelling that the polynomials are now restricted to only the four variables 
x, I, z,n and we can put the polynomials p 2 = q% — 0. This simplifies the equations §3§ 
considerably. The second and third of the equations in ^ simplify to: 

92 = -Pll, 

P3 = -qiz. (6) 

Substituting these in the first equation, we need to solve just one equation: 

Pl (l - Ix) - qi (l - zn) = , (7) 

to obtain the full polynomial vector as a solution. It is clear that solutions are of the Michelson 
type denoted by X. One such solution has already been given in the literature [21 [B]. This 
solution in our notation is: 

Pi = 1 — zn — znlx + lx(zn) 2 

qi = 1 — Ix — Ixzn + zn(lx) 2 . (8) 
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Writing, 

A = Pl (l - Ix) - gi {l - zn) , (9) 

we get for (|5J), A = [znlx,lxzn] which is a commutator that vanishes. Thus it is an element of 
U and ([HI) is a solution (over the quotient ring). This is already well known. 

What we would like to emphasise is that there are more solutions of this type - in fact there 
are infinite number of such solutions. In what follows we give a systematic way to enumerate 
and list such solutions. 

Notice that the operators Ix and zn occur together. Physically, these operators correspond 
to round trips. We denote these round trip operators by a = Ix and b — zn, then we can write 
the previous equation and solution respectively, in a compact way as: 

pi(l-o)-«i(l-6)=0, (10) 

and, 

Pi = 1 — b — ba + ab 2 , 

qi = l- a-ab + ba 2 . (11) 

We can substitute this in A = P i(l — a) — qi(l — b) to obtain A = [ba,ab] which vanishes and 
hence (fTTj) is a solution. 

But we can proceed further and write: 

Pi = Pi~ ba 2 b , 

q[ = qi - ah 2 a , (12) 

to obtain: 

A' = ba 2 ba - ab 2 ab. (13) 

A' however does not vanish - it is also not a commutator - and hence p[ , q[ is not a solution. 
Continuing further with, 

Pi — Pi + ab 2 ab , 

q'l = q[+ ba 2 ba , (14) 

we get, A" = ba 2 bab — ab 2 aba which also does not vanish, again not a commutator. 

This process can be continued in such a way so that lower degree terms cancel and so as to 
get another vanishing commutator but of higher degree. The next commutator is obtained by 
adding two more terms of successively higher degree in a similar way. The result is: 

pi = 1 — b — ba + ab 2 — ba 2 b + ab 2 ab + ab 2 aba — ba 2 bab 2 , 

qi = 1 — a — ab + ba 2 — ab 2 a + ba 2 ba + ba 2 bab — ab 2 aba 2 . (15) 

Note that qi is obtained from p\ by interchanging a and b and vice versa. This is a symmetry. 
The expressions for p\ and qi in Eq. (|15p give the vanishing commutator: 

A = [ab 2 a, ba 2 b] , (16) 

and hence (fl"5"j) is a solution. 

There are two more solutions with eight terms. One of them is: 

pi = 1 + a - b 2 - b 2 a - b 2 a 2 - b 2 a 3 + a 2 b 4 + a 2 b 4 a 

= (1-b 2 -b 2 a 2 +a 2 b i ){l + a), 
qi = 1 + b - a 2 - a 2 b - a 2 b 2 - a 2 b 3 + b 2 a 4 + b 2 a 4 b 

= {l-a 2 -a 2 b 2 + b 2 a 4 )(l + b), (17) 
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with A = [b 2 a 2 ,a 2 b 2 ] which is also a vanishing commutator. This solution looks something like 
the square of the solution obtained by Tinto et al (see [5]) multiplied by linear factors in 
order ensure successive term cancellation. Here one starts with p\ — 1 + a... and continues until 
one reaches the above commutator. 

The third solution corresponds to the commutator [baba, abab]. This solution is given by: 



The following observations are in order: the solutions correspond to the vanishing 
commutators as follows; there is only one commutator at degree 4, namely, [ba,ab]; of degree 
8, there are three commutators: [a 2 b 2 ,b 2 a 2 ], [abab, baba], [ab 2 a, ba 2 b}. And these are the only 
three commutators at degree 8. Also there are equal number of a's and b's in each string of the 
commutators. To each of the vanishing commutators, there is a corresponding polynomial vector 
solution (pi,<7i). The rest of the solution vector is obtained from Eq.© which then yields p% 
and q 2 . 

These are solutions only upto degree 7. Clearly one can construct solutions of higher degree. 
We do this in a combinatorial way in the next section. We enumerate the solutions in a given 
order - called the lexicographic order - and also present an algorithm to obtain the polynomials 
Pi and q\. 

4. Solutions of higher degree 

4-1- Listing the commutators 

In order to proceed further, let us first define an interchange operator / between the symbols a 
and b, that is, given a string s composed of a and b, the string t = Is is obtained by replacing a by 
bin s. A string is a polynomial with only one term usually also called a monomial in the literature. 
The length of a string is the degree of the monomial. We will use these terms interchangably, 
but that should not lead to any confusion. Note that all the vanishing commutators which have 
appeared corresponding to a solution have the form [s,Is]. Also since the string Is must be a 
permutation of the string s, in order that the commutator [s,/s] vanish, there must be an equal 
number of a's and 6's in the string s of the vanishing commutator. So if we have n number of a's 
in the string s, it must also have the same n number of b's, so that the total string s has length 
2n. The commutator is then made up of strings of length An. Note that the operator / can be 
extended to a polynomial containing strings in an obvious way by defining: 



Then we easily check that in the foregoing qy — Ip\, where (jpx,qi) constitute a solution to 
Eq. ([TP)). Also I 2 = 1 the identity mapping. / can be considered to be an operator over the 
polynomial ring over the indeterminatcs a and b. 

In the foregoing we have then considered cases for n = 1,2. For n = 1 there is one solution, 
while for n = 2 there are three distinct solutions. In order to proceed further, we must first list 
the commutators at each n. In order to carry out this programme, we must introduce an order 
among the strings of length 2n. We choose this order to be lexicographic or if we consider strings 
of all lengths, then length lexicographic [S]- But since at a time we deal only with strings of a 
given length 2n, the lexicographic order suffices. We arbitrarily choose a < b. Then at length 
2, aa < ab < ba < bb. At length 3, aaa < aab < aba < abb < baa < bab < bba < bbb and so 
on. More generally, for strings of the same length, s < t if at the first position from the left the 
strings differ, the symbol in string s is less than the symbol in string t. In this way we may list 
all strings in the symbols a, 6 of a given length. This is all we require. 



pi = 1 — b + ab — bab — baba + abab 2 — baba 2 b + abab 2 ab , 
qi = 1 — a + ba — aba — abab + baba 2 — abab 2 a + baba 2 ba . 



(18) 



I(s 1 + s 2 + ...) = I Sl +Is 2 + .... 



(19) 
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Thus at length 2n = 2 we have just one commutator [ab, ba]. At 2n = 4, we have 
[aabb,bbaa],[abab,baba],[abba,baab], that is, three commutators; we have also three solutions 
corresponding to each commutator. We have also taken care to write the first string in the 
commutator to be less than the second string in the lexicigraphic order chosen. Note that we 
could list the commutator just by the first string, because the second string in the commutator 
is obtained by applying the operator / to the first string. We will give an algorithm to obtain a 
solution, that is, the polynomials p\ and gi in a unique way given the commutator. Thus listing 
all the commutators of a string of length 2n amounts to listing the solutions (pi,qi) having that 
commutator. The polynomials pi and qi have 2n terms and degree 2n — 1. We are now ready to 
list the commutators [s,Is] with first string s of length 2n. 

Consider n = 3. The string s has length 6 with 3 a's and 3 o's. Listing the strings in 
lexicographical order we have: aaabbb < aababb < aabbab < aabbba < abaabb < ababab < 
ababba < abbaab < abbaba < abbbaa. We stop here because, the next string is baaabb = I{abbbaa) 
and abbbaa has occurred in the listing before and hence would lead to minus the commutator 
and therefore not an independent commutator. This follows from the symmetry encoded in the 
interchange operator /. Thus at n = 3 there are 10 independent commutators. Proceeding in this 
way we find that at n = 4 there are 35 such commutators. The general formula for the number of 
commutators is C n -\. This formula follows from the following argument. Consider a string 
of length 2n, with n number of a's and 6's. Because of symmetry between a and b, the first a is 
fixed in the string; as seen, the string starting with b does not lead to a new commutator. Thus 
we must choose n — 1 combinations of a's among 2n — 1 available positions in the remaining 
string. This is the number of independent commutators which is then just n C n -\. 

4-2. Algorithm for the solution given the commutator 

The algorithm is obtained by observing the way we got the solutions of degree 7 for n = 2. We 
started with the lowest degree and proceeded towards obtaining the commutator by cancelling 
terms at successively increasing degree. But now we have been assigned a commutator before 
hand and we need to find the solutions (pi,qi). We therefore need to go in the reverse order. 
We freely make use of the operator / in this algorithm. 

Consider the commutator A = [s 2n ,Is2n\ = S2 n I{s2n) — I{s2n)s2n where s 2n is a string 
containing n number of a's and an equal number of 6's. 

(i) Consider the first string in the commutator A which we denote by t^ n — S2 n I(s2n)- The 
subscript denotes the length of the string or the degree of the monomial. Now, either t 4n 
ends in a or b; that is, t 4n — Si n -ia or t 4n = Si n -\b. 

• If t 4n — S4 rl -ia, then write i4„_i = S4„-i or; 

• if U n = Sin-ib, then write i 4 „_i = -7s 4 „_i. 

(ii) Replace ti n by ti n -i and repeat. That is either we have ti n -\ = Si n -2<i or ti n -\ = S4n-20- 
If U 

n—i — s in—2 a : then write ti n —2 — s in— 2 otherwise write ti n —2 — Isin—2- 

(iii) Continue this process till one arrives at a string of length 0, namely, t . Note that t = ±1- 

(iv) Then p 1 = Y^t^ and °i = I P^- 

Note that tk are strings of length k. p\ and q\ are polynomials of degree An — 1. 

To illustrate the algorithm explicitly, let us take the case for n = 3 and choose one 
commutator among the 10 possible ones. Let us randomly choose the string ababba, so that 
the commutator is [ababba, babaab}. So we have sq — ababba. Thus t 12 = abab 2 ababa 2 b. Now 
t\2 ends in b. Thus we write t\ 2 = sub, where sn = abab 2 ababa 2 . Thus choosing the second 
option in step (ii), we have in = —Isu = —baba 2 babab 2 . This is the highest degree term in the 
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polynomial p\ and for n = 3, it is of degree 11. The next term is obtained from tn. Again t\\ 
ends in b. Following the steps of the algorithm, tn = siob and hence tio = —Isiq = abab 2 ababa. 
Since, tio ends in a, tg = sg = abab 2 abab. We then have t% — —baba 2 bab and so on. We must 
carry out the steps till we reach degree 0. This terminates the algorithm. In this way, we can 
easily construct the 11th degree polynomial p\ by adding up all the strings tk, k = 1, 2, 11. 
Collecting all the terms and writing out the polynomial p\ explicitly: 

pi = 1 — b + ab — bob — baba + abab 2 — baba 2 b + abab 2 ab 

— baba 2 bab + abab 2 abab + abab 2 ababa — baba 2 babab 2 . (20) 

The polynomial q\ is obtained from p\ by just interchanging a and 6, that is, q\ = lp\. The 
other polynomials q2 and p% are obtained from Eq. ((5]). These polynomials then constitute the 
polynomial vector which is a solution and is in the module of solutions over the quotient ring 
K. In terms of the time-delay operators x, I, z, n, the polynomials p\ and q\ are of degree 22, 
while q2 and pg, are polynomials of degree 23. For a general value of n, the solution contains 
polynomials of maximum degree 8n — 1 in the time-delay operators. 

Clearly, from the mathematical point of view we have an infinite family of solutions. Note 
that here we do not lay any claim to any exhaustive listing of solutions. However, we believe this 
family of solutions is sufficiently rich, because we can form linear combinations of these solutions 
and then these form a left submodule of the full left-module of solutions. We can also subtract 
solutions like (jlip from (|15p which is also a solution with each component polynomial having 
less number of terms (or links in the usual terminology) albeit of the same degree. Further, by 
using cyclic permutations on {x, y, z) and (/, m, n) we can generate more solutions. For example, 
we may set a — my and b — zn and obtain solutions for the four links, y,z,m,n. The linear 
combinations of these solutions with the previous ones in the variables x, z, I, n are also solutions 
in the full module pertaining to all six links. Thus we see that a plethora of solutions can be so 
generated. 

From the physical point of view, since we have neglected terms in L and L 2 and higher 
orders in our calculations, a limit on the degree of the polynomial solutions arises. That is upto 
certain degree of the polynomials, we can safely assume the commutators to vanish. But as the 
degree of the polynomials increases it is not possible to neglect these higher order terms and 
then such a limit becomes important. We now investigate this limit and make a very rough 
estimate of it. Our LISA model in (8] gives L ~ 10~ 6 m/sec 2 . The other term L 2 /c 2 ~ 10~ 15 is 
way too small to matter. So the limit is set by L. From L we compute the error in L, namely, 
AL ~ ^At 2 L. If we allow the error to be no more than say 10 metres, then we find At ~ 4500 
sec. Each time-delay is about 16.7 sees which gives the number of successive time-delays to be 
about 270. This is the maximum degree of the polynomials. This means one can go upto n < 30. 
If we set the limit more stringently at AL ~ 1 metre, then the highest degree of the polynomial 
reduces to about 80 which means one can go upto n = 10. This means there are a large number 
of TDI observables available to do the physics. We however point out that when we go to higher 
degree polynomials, the number of interpolations of data also increases leading to accumulation 
of numerical errors. This is a relevant issue whenever high degree polynomials are used. 

5. Conclusion 

In this work we have given a family of approximate solutions for the general model of LISA where 
the up-down links are unequal and the armlengths are time varying, but with only two arms 
functional. Such a situation can arise when one of the arms is non-functional or when it is not 
possible to obtain data from all optical links because high costs etc. We then have only four links 
to deal with and although the general non-commutative problem of obtaining second generation 
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TDI is extremely difficult, the problem here reduces to a much simpler and tractable one. We 
have therefore been able to exhibit a family of solutions although not exhaustive, is quite rich 
and interesting both from the mathematical and physical point of view. Also the mathematical 
analysis presented here could be extended to other space missions such as ASTROD where a 
similar situation arises although the physical parameters will be quite different. 
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